Critical parameters for the disorder-induced metal-insulator transition in FCC and 

BCC lattices 
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We use a transfer-matrix method to study the disorder-induced metal-insulator transition. We 
take isotropic nearest- neighbor hopping and an onsite potential with uniformly distributed disorder. 
Following previous work done on the simple cubic lattice, we perform numerical calculations for the 
body centered cubic and face centered cubic lattices, which are more common in nature. We obtain 
the localization length from calculated Lyapunov exponents for different system sizes. This data 
is analyzed using finite-size scaling to find the critical parameters. We create an energy-disorder 
phase diagram for both lattice types, noting that it is symmetric about the band center for the 
body centered cubic lattice, but not for the face centered cubic lattice. We find a critical exponent 
of approximately 1.5-1.6 for both lattice types for transitions occurring either at fixed energy or at 
fixed disorder, agreeing with results previously obtained for other systems belonging to the same 
orthogonal universality class. We notice an increase in critical disorder with the number of nearest 
neighbors, which agrees with intuition. 



I. INTRODUCTION 

The disorder-induced metal-insulator transition (MIT) 
and the concept of Anderson localizationii^ii^ii^i^ have 
been studied extensively for more than forty years. The 
scaling theory of localization® provides a very success- 
ful approach for non-interacting electrons. According to 
its predictions, a disorder driven MIT occurs in three- 
dimensional systems, i.e. beyond a critical amount of dis- 
order Wc all eigenstates localize. For smaller disorder, 
extended states exist in the system. 

For the simple cubic (SC) lattice and uniform dis- 
order distribution, the critical disorder and the crit- 
ical exponent have been successfully calculated using 
the transfer- matrix method {TMM)Jj^ Highly accu- 
rate recent studies^ii^iiiii^ report Wc = 16.54 ± 0.02 
and the critical exponent v — 1.57 ± 0.02.^ How- 
ever, direct diagonalization results based on energy-level 



statistics 
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and multi-fractal analysis 



,17,18 



give a 



smaller v = 1.44±0.2. Furthermore, experimental results 
report yet smaller values oi v > 1 .0 ,il22j^L22j21 Many of 
these observed discrepancies can be explained by the at- 
tainable limits on system sizes, temperature and statisti- 
cal averages in the above results, nevertheless, the quest 
for an accurate determination of the critical parameters 
at the Anderson transition is not yet complete. 

The SC lattice is, in addition, not very common?4, in 
nature. The only element known to adopt it is the alpha 
phase of polonium; most metals exhibit body centered 
(BCC) or face centered (FCC) cubic lattices. Although 
the localization properties for an FCC lattice have been 
studied recently^^ for a vibrational problem, to the best 
of our knowledge no critical parameters have been re- 
ported for an electronic Anderson transition in BCC and 
FCC lattices. In order to fill this gap, we use in the 
present paper the TMM and finite-size scaling (FSS) to 



calculate the critical parameters for the MIT in the BCC 
and FCC lattices. Of course, one should expect no change 
in the critical exponent, as all systems considered here 
belong to the same orthogonal universality class. On the 
other hand, because of the different numbers of nearest 
neighbors in SC {Z = 6), BCC (8) and FCC (12) lattices, 
the values of the critical disorder Wc and energy Ec may 
be different. Hence our study tests and reconfirms uni- 
versality while at the same time allowing to see how the 
non-universal parameters of the transition change with 
increasing coordination number. 



II. NUMERICAL APPROACH 

A. The Transfer-Matrix Approach to 
the Anderson Model of Localization 

To model the MIT in the 3D system we use the stan- 
dard Anderson Haniiltonian 



(1) 



The orthonormal states \i) correspond to electrons lo- 
cated at sites i = [x, y, z) of a cubic lattice with pe- 
riodic boundary conditions. The hopping integrals tij 
are non-zero only for i, j being nearest-neighbors and the 
energy scale is set by choosing tij = 1. The disorder 
in the model is incorporated into the diagonal energies 
Ci S [—W/2, W/2], randomly distributed according to the 
uniform distribution with width W . 

In order to compute the localization length A of the 
wave function, we use the TMM for quasi-lD bars of 
cross section M x M and length L > MJ^i^^^ The 
Schrodinger equation Hip = Eijj for the Hamiltonian 



given by Eq. ([T]) is written in the TMM form: 



-I 



1 



HO 







ipi 



where "^^H; and T; denote the wave function, Ha: 
tonian matrix and transfer matrix of the /th shoe of 
bar, respectively. 1 and denote unit and zero m; 
ces. The locahzation length X{M,W) = l/7min at 
ergy E is determined by the smallest Lyapunov e? 
nent 7min > obtained as an eigenvalue of the proc 
of transfer matrices tl = TiTi_i...T2Ti, where L is 
increased until the desired accuracy is achievedi^^ The 
reduced localization length may then be calculated as 
Am{W) = \{M, W)/M. 

Ci and Ci_|_i are the connectivity matrices describing 
the connections of the Zth slice to slices I — 1 and / + 1.^^ 
Element Cjk of the connectivity matrix equals 1 if the 
site j in one slice is connected to the site k in the other; 
otherwise Cjk = 0. In the case of the SC lattice, each site 
has only one connection to the succeeding (preceding) 
layer; therefore all Ci are unit matrices and the transfer 
matrix T; reduces to the most often used forro^^ 
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For BCC and FCC lattices the connectivity matrices take 
a more complicated form, but with purely diagonal disor- 
der, i.e. no disorder in the hopping integrals Uj, they are 
constant so that the inverse Cf^ needs to be calculated 
only once at the beginning of the TMM calculations for a 
given size M. Nevertheless, the additional need to multi- 
ply all states at each step of the TMM with a dense ma- 
trix reduces the speed of the calculation and hence 
restricts the attainable system sizes. 



B. The Lattice Structures 

The structure of the BCC lattice is displayed in Fig. 
[TJa). The construction of the TMM quasi-lD bar pro- 
ceeds along a (100) vector. In this case each site within 
the slice is connected to four sites in the preceding slice 
and to four sites in the succeeding one. There are no 
connections between sites within the slice, which means 
that the Hamiltonian matrix H; is a diagonal matrix of 
energies e^. We use periodic boundary conditions in both 
transversal directions, which results in the connectivity 
matrix for a slice of M x M sites being singular for all 
even M, thus restricting the system sizes we can use. 
Using a helical boundary condition^— in one or two direc- 
tion provides the same singularities and hence offers no 
advantage. 

Fig. [TJb) shows the structure of the FCC lattice. It 
proved convenient to construct the TMM bar along a 




FIG. 1: (a) 3 layers of the 3D BCC lattice along a (100) 
lattice vector. Light gray spheres mark the 1st and 3rd layer, 
dark gray ones indicate the central layer. Lines between layers 
denote the connections between lattice sites. The connections 
to the upper-left sphere in the central layer are emphasized 
by broad lines, illustrating its 8 neighbors. The 4 thick light 
gray lines connect 1st and the central layer, the black ones go 
from the central to the 3rd layer, (b) Structure of 3 layers of 
the 3D FCC lattice along a (111) lattice vector. The broken 
lines mark the cubic unit cell of the lattice. Sites in the 1st 
and 3rd layer are dark gray; sites in the central layer are 
light gray. The thick light gray lines represent connections 
between lattice sites in the same layer for one particular site. 
The thick black lines represent connections between lattice 
sites in neighboring layers for another site. The thin lines 
indicate connections to other sites. Some sites in the upper- 
right corner are removed for clarity. 



(Ill) vector, so the subsequent layers of the bar are close 
packed. Within the layer each site has six connections 
to nearest neighbors. In addition there are three con- 
nections to the preceding and three connections to the 
succeeding layer. The resulting connectivity matrix can 
be inverted for each size of the M x M TMM slice but 
only when we use a mix of periodic boundary conditions 
in one direction and helical boundary conditions in the 
other. See appendix for examples of the connectivity ma- 
trices for system size M = 3. 



C. Finite-Size Scaling 

The MIT is characterized by a divergent correlation 
length, so that at fixed energy E, ^{W) oc \W - Wd'" 
and at fixed disorder W, ^{E) oc - Ec\~'^, where v is 
the critical exponent and Wc, Ec are the critical disorder 
and energy, respectively, at which the MIT occurs.^ In 
the following discussion we shall assume the case of fixed 
energy and varying disorder; the converse case of fixed 
disorder and varying energy proceeds analogously. 

In order to extract the critical parameters from the 
calculated values of Am( W), one applies the FSS proce- 
dure outlined in Ref. [30. The correlation length for the 
infinite system ^ may be obtained from the localization 
lengths for finite system sizes Am{W) by using the one- 
parameter scaling law Am = f{M/(^)M- The FSS can 
be performed numerically by minimizing the deviations 
of the data from a common scaling curve. The criti- 
cal parameters are then obtained by fitting the ^ values 



as obtained from FSS. Better numerical accuracy for the 
FSS procedure can be achieved by fitting directly the raw 
data from TMM calculations using the method applied 
previously to the TMM data for the 3D SC lattice,^ii^ We 
introduce a set of fit functions which include two kinds 
of corrections to scaling, (i) nonlincarities of the W de- 
pendence of the scaling variables and (ii) an irrelevant 
scaling variable which accounts for a shift of the point at 
which the Am{W) curves cross. We usei^ 



AM-/(XrM'/^XiM^) 



(4) 



where Xr and Xi the relevant and irrelevant scaling 
variables respectively. The function Am{W) is then Tay- 
lor expanded 



Am = Y.x^M^yUxrM'/-') , (5) 

rir 

fn = • (6) 



k=0 



Nonlincarities are taken into account by expanding Xi 
and Xi in terms oi w — {Wc — W)/Wc up to order mr 
and TTii, respectively, 
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m=0 



with &i = Co = 1 . The expansions in the fit functions are 
carried out up to orders ni, rir, Wp, m-i which are adjusted 
to the specific data and should be kept as low as possible, 
while giving the best fit to the data and FSS plot and 
minimizing the errors for critical parameters and v. 
The Levenberg-Marquardt method was used to perform 
the non- linear fitJ^i^^. 

We emphasize that this FSS procedure assures the di- 
vergence of ^ and hence it is not the divergence itself but 
rather the quality of how the model fits the computed 
reduced localization lengths Am which determines the 
validity of the scaling hypothesis. 



III. CALCULATIONS AND RESULTS 
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FIG. 2: Phase diagram for the BCC lattice. The dark gray re- 
gion represents the approximate location of the phase bound- 
ary. Its edges (the solid black lines) were determined by com- 
paring localization lengths with errors < 10% for system sizes 
M = 7 and M = 9 in the [E, W) plane. The solid squares 
(□) are points calculated by performing high-precision FSS on 
localization data with an error < 0.1%. The dashed squares 
are reflections of the solid squares in the i? = axis. The 
diamonds (♦) denote the band edges at = 0. They have 
been joined to the phase boundary edges calculated for higher 
disorders as a guide to the eye. The dashed lines are the the- 
oretical band edges ±{Z -\- W/2), where Z is the coordination 
number. The light gray, shaded area in the center contains 
extended states; states outside the phase boundary are local- 



A. Phase Diagrams 




Fig. [2] and Fig. [3] show the phase diagrams for the 
BCC and FCC lattices, respectively. Originally a grid 
of W versus E values was created with separation AE, 
AW = 0.5. At each point the nature of the elec- 
tronic wavefunction was determined by comparing the 
reduced localization lengths Am calculated for system 
sizes M = 7 and M = 9 with error < 10%. If Ag > Ay 
(<) at the same values of E and W then we identify the 
point {E, W) in the phase diagram as extended (local- 
ized). The edges of the phase boundary were obtained 
by averaging separately over the three extended and lo- 
calized points {E, W) nearest to the boundary and then 



FIG. 3: Phase diagram for the FCC lattice. Symbols, lines 
and shaded areas have the same meaning as in Fig. [5] with the 
diamonds (♦) representing the band edges —12 and 4 at zero 
disordes^i. Error bars are within symbol size. 



connecting such averages using a spline fit. We do not 
obtain data points for lower disorder values, as the fiuc- 
tuations in the Lyapunov exponents, due to the small 
system sizes, become too big; higher values of disorder 
smooth out these fluctuations. 

A striking difference between the phase diagrams is 
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FIG. 4: Reduced localization lengths Am versus disorder W 
for BCC lattice. System sizes M are 3(«), 5(0), . . . , 15(x). 
Error bars are within symbol size. Lines are fits to the data 
given by Eqs. (0} - ((Tjl with rir = 3, rii = 2, mr = 3, mi = 1. 



that for the BCC lattice the phase boundary is symmet- 
ric about the hnc E = 0, whereas for the FCC lattice it 
is not. This is due to the bipartiteness of the BCC lat- 
tice which consists of two SC sublattices, one displaced 
half the distance along a body diagonal of the other. 
Hence for any site in one sublattice, its nearest neigh- 
bors are in the other sublattice. Such connections result 
in states coupled by a bipartite symmetry transformation 
— which is exact for the case of no diagonal disorder — 
with eigenenergies of the same magnitude but opposite 
sign having approximately the same localization lengths; 
this produces a symmetric phase diagram. The FCC lat- 
tice is non-bipartite, so such a symmetry in its phase 
diagram is not observed. 



B. Critical Parameters at i5 = 

The TMM calculations were performed for system sizes 
up to M = 15. In order to examine the localization 
properties at the band center for the BCC lattice and 
the barycenter— for the FCC lattice, we set £' = in 
Eq. A value of the critical disorder Wc was approx- 
imated using the phase diagrams described above and 
then localization lengths A were calculated for a range of 
W close to this approximate value with accuracy ranging 
from 0.1% for small system sizes M to about 0.14% for 
the largest. Let us remark that we use the term critical 
disorder to indicate that there are no further extended 
states at E' = for disorders W > Wc', extended states 
may still exist for W > Wc at other energies E, as shown 
in Fig. O 

The reduced localization lengths for the BCC lattice 
are displayed in Fig. [4l Note how the crossing point of 
the curves shifts with changing AI. In most cases this 
indicates the need for an irrelevant scaling variable in- 
troduced via non zero values of rii and mj in Eqs. ([5]) 
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FIG. 5: Scaling function (solid line) and scaled data points for 
the BCC lattice and rir = 3, ni = 2, nir = 3, mi = 1. Symbols 
denote the same values of M as in Fig. |4] Inset: Dependence 
of the scaling parameter ^ on the disorder strength W for the 
1 Q ix^ ,,„i..„„ r?;„ [71 T„ „n „„„„„ „.;+u;„ 
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FIG. 6: Reduced localization lengths Am versus disorder W 
for FCC lattice. Symbol sizes M are 3(»), 4(v), . . . , 15(x). 
Error bars are within symbol size. Lines are fits to the data 
given by Eqs. (H} - ((TJ with rir = 2, rrir — 2. Lines for even 
M have been removed for clarity. 



and ([7])- Fig. [5] shows the results of the scaling proce- 
dure for Tij- = 3, Tij = 2,TOr = 3, mi = 1. The scaling 
curve exhibits localized and extended branches as ex- 
pected for the MIT. Divergence of the scaling parameter 
^ a.t W 20.75 indicates the critical value of the disor- 
der. Table llja) gives some examples of models providing 
the best fits and the resulting critical parameters. The 
values of the critical disorder and critical exponent ob- 
tained by averaging over all the best fit models are also 
given. 

Results of the TMM calculations for the FCC lattice 
are shown in Fig. [S] In this case the lines for constant M 
cross at the same point — at least within the accuracy 
of the calculated Am — indicating that the use of the 



TABLE I: Critical parameters for the MIT in the BCC lat- 
tice. All errors quoted are standard errors, (a) 3 examples of 
FSS results with varying Tir, ni, mr, mi at fixed energy E — 0. 
We use 91 data points, equally spaced in the indicated in- 
tervals (cp. Fig. |4}, for each set of nr,ni,mr,mi. Varying 
nr, rii, mr, mi, we obtain 77 best fit models in order to pro- 
duce the indicated averages, (b) Similar FSS results obtained 
for 3 out of 19 best fit models from 82 non-equally spaced data 
points at fixed W = 15 for the indicated energy intervals, (c) 
Results at W = 17.5 (cp. Fig. [8]) for 3 out of 8 best fit mod- 
els with 108 non-equally spaced data points used in each FSS 
procedure. The numerical fitting procedure continued in all 
cases until convergence was reached or (a) 5000, (b,c) 1000 
iterations had been completed. 

(a) 



AM E AW Wr m rur rm Wc u y_ 

3 - 15 20.3 - 21.5 2 1 20.95(1) 1.67(5) 

3 - 15 20.3 - 21.5 3 1 1 4 20.92(2) 1.51(9) 1.7(5) 

3 - 15 20.3 - 21.5 3 2 3 1 20.75(3) 1.70(9) 3.0(5) 



averages: 20.85(1) 1.61(1) 



(b) 



AM AE W Ur m rm Ec ^ y 

9 - 13 9.9 - 10.9 15 2 1 10.38(1) 1.32(5) - 

9 - 13 9.9 - 10.9 15 2 2 10.38(1) 1.22(5) - 

9 - 13 9.9 - 10.9 15 3 4 10.40(1) 1.03(3) - 



averages: 10.39(1) 1.27(1) 

(c) 



AM AE W rir rn mi Ec i' y 

7 - 15 10.5 - 11.5 17.5 2 1 10.98(1) 1.55(6) - 

7 - 15 10.5 - 11.5 17.5 3 2 10.99(1) 1.48(6) - 

7 - 15 10.5 - 11.5 17.5 3 4 10.99(1) 1.36(7) - 



averages: 10.96(3) 1.47(9) 



irrelevant variables in Eqs. ^ and ([7]) is not necessary 
in most cases and n-i = m-i ~ Q. Results of the fit for 
rir = 2, TOi- = 2 are displayed in Fig. [71 The transition at 
W ~ 26.73 is clearly indicated. More examples of best fit 
models can be found in Table|TT](a) as well as the average 
values of the critical parameters. 



C. Critical Parameters away from tiie Band Centre 

We also perform calculations where we fix the disorder 
and allow the energy to vary across a critical value for 
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FIG. 7: Scaling function (solid line) and scaled data points 
for the FCC lattice and rir = 2, rrir = 2. Symbols denote 
the same values of M as in Fig. [6] Inset: Dependence of the 
scaling parameter ^ on the disorder strength W for the 13 
W values shown in Fig. [6] In all cases error bars are within 
symbol size. 



TABLE II: Critical parameters for the MIT in the FCC lat- 
tice. The fitting procedure was continued until convergence 
was reached or until 1000 iterations had been completed, (a) 
At _E = (cp. Fig. O, we use 91 data points for each FSS 
and the obtained 31 best fit models average as indicated, (b) 
83 data points (cp. Fig. [9]) are used to perform the FSS at 
= 18 and 8 best fit models are used to obtain the aver- 
ages. 

(a) 



AM E AW Ur Hi mr rm Wc v y 

3 - 15 26 - 27.5 1 4 3 4 26.72(1) 1.49(4) 6(2) 

3 - 15 26 - 27.5 2 2 26.73(1) 1.58(3) 

3 - 15 26 - 27.5 3 4 2 4 26.72(1) 1.51(10) 4(2) 

averages: 26.72(1) 1.53(1) 

(b) 

AM AE W m ni m^ m-i Ec v y 

9 - 15 8.52 - 8.88 18 1 2 8.683(3) 1.63(5) - 

9 - 15 8.52 - 8.88 18 2 2 8.687(4) 1.65(5) - 

9 - 15 8.52 - 8.88 18 3 1 8.685(3) 1.62(5) - 



averages: 8.684(2) 1.63(2) 



the transition. We remark that it is know that such in- 
vestigations are numerically more difficult due to the in- 
fluence of density of states effects." Results of the TMM 
and FSS calculations for the BCC lattice with W = 11. b 
can be seen in Fig. [5] Although the lines for constant 
M do not all cross at the same point, the data is best 
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FIG. 8: Localization data for the BCC lattice with W = 17.5. 
System sizes M are 7(a), 9(+), . . . , 15(x) as in Fig. g] Error 
bars are within symbol size. Left: Scaling function (solid line) 
and scaled data points using rir = 2,mr — 1. Right: Reduced 
localization lengths Am versus disorder W. Lines are fits to 
the data given by Eqs. © - (O with = 2, mr = 1. 
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FIG. 9: Data for the FCC lattice with W = 18. System 
sizes M are 9(+), 11(0), • • • , 15(x) as in Fig. [7] Error bars 
are within symbol size. Left: Scaling function (solid line) and 
scaled data points using rir — l,mr — 2. Right: Reduced 
localization lengths Am versus disorder W . Lines are fits to 
the data given by Eqs. - ([7]) with rii — 1, rrir = 2. 



fitted by models containing no irrelevant scaling. Also 
evident is the poorer quality of the fit compared to the 
calculations where the energy was fixed at zero. This is 
not due to using lower accuracy data, as the maximum 
raw-data error remained at 0.1%. Hence we attribute it 
to complications arising from a varying density of states 
close to Ec at the attainable values of M. Results for the 
critical parameters are shown in Table [J (b) and (c) for 
W — 15 and W — 17.5 respectively. The low value of 
for the case W — 15 can be attributed to the use of fewer 
data points in the FSS and only using three values of M. 
We note that this is consistent with the lower values of 
v obtained in the diagonalization studies as mentioned 
in the Introduction. It appears that the FSS procedure 
systematically reduces the values of the critical exponent 
for data from smaller systems or of lower accuracy. 

Resuhs for the TMM and FSS calculations for the FCC 
lattice with W = 18 can be seen in Fig. M Table [IT] (b) 
gives examples of the best fit models and shows the re- 
sulting average critical parameters. Note that both esti- 
mates of v are consistent with the result 1.57(2) for the 
SC lattice,^ 



IV. CONCLUSIONS 

Using the transfer-matrix approach and FSS we deter- 
mined the critical parameters of the Anderson transition 
for the BCC and FCC lattices. The values of the crit- 
ical exponent v are in good agreement with the results 
obtained previously for other systems belonging to the 
orthogonal universality class. The increase of the crit- 
ical disorder Wc from 16.54 for the SC lattice to 20.85 
for BCC and 26.72 for FCC lattice may be attributed 
to an increasing number of nearest neighbors which for 
the above structures equals 6, 8 and 12, respectively. 



More nearest neighbors connected to a given site provide 
more paths for electronic transport, so stronger disor- 
der is needed to localize eigenstates of the system. The 
universal localization properties of a 3D system and the 
presence of an MIT are however not affected in accor- 
dance to the scaling theory of localizatior^ and in agree- 
ment with results21 showing that they depend only on 
the dimensionality of the system, but not on the number 
of nearest neighbors in the lattice. 

Our results and their interpretation are consistent with 
investigations of classical bond and site percolation mod- 
els on SC, BCC and FCC lattices. In Ref . Issl it was found 
that the percolation thresholds for these lattices decrease 
with increasing number of nearest neighbors; more neigh- 
bors allow for easier formation of a percolating cluster, 
or, as in our case, the formation of extended states. 
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APPENDIX A: CONNECTIVITY MATRICES 

For completeness, let us give the connectivity matrices 
for BCC and FCC lattices with M = 3. Recall that 
C/ is the connectivity matrix describing the connections 
of the Ith slice to the I — 1th slice. Element Cjk of the 
connectivity matrix equals 1 if site j in the ^th slice is 



7 



connected to site k in the I — 1th sHce; otherwise Cjk = 0. 
The boundary terms are indicated in itaUcs. For the BCC 
lattice for odd layers, 



For the FCC lattice for odd and even layers 



For even layers 



/I 


1 











1 





1 \ 


























1 1 














1 


1 





























1 


1 














1 


1 








/ 1 


1 

















1 


0\ 


1 


1 


1 





1 




















1 


1 

















1 


1 1 





1 


1 
















(Al) 




1 





1 


1 

















1 


1 





1 


1 




















1 





1 


1 




















1 





1 


1 





1 






C/ = 








1 





1 


1 

















1 


1 





1 


1 




















1 





1 


1 








\0 








1 


1 





1 


I J 




















1 





1 


1 










































1 





1 


1 
























\1 

















1 





1/ 



(A3) 



/I 


1 





1 


1 o\ 





1 


1 





110 


1 





1 


1 


10 











1 


10 110 














110 11 











1 


1 i 1 


1 


1 










110 





1 


1 





1 1 



(A2) 



In all cases, I = 1,2, 
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